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1 INTRODUCTION 



ABSTRACT 

We present extended simulations of the large-scale distribution of galaxy clus- 
ters in several different dark matter models, using an optimized version of the 
Zel'dovich approximation. The high accuracy and low computational cost of this 
technique allow us to run a large ensemble of fifty realizations for each model, 
and we are therefore able to make an accurate determination of the cosmic 
variance. Six different dark matter models are studied in this work: Standard 
CDM, Open CDM (fio = 0.2), Tilted CDM with primordial spectral index 
n = 0.7, mixed Cold + Hot dark matter with fJhot = 0.3, spatially flat low- 
density CDM with cosmological constant term corresponding to Q^ = 0.8 and 
low Hubble constant [h = 0.3) CDM. We compare the discrete and smoothed 
cluster distributions with those of an Abell/ACO redshift sample using various 
statistical quantities, such as the J3 integral and the probability density func- 
tion (pdf). We find that the models that best reproduce the behaviour of J3[R) 
(i.e. the integral of the two-point correlation function), the pdf and the mo- 
ments of the Abell/ACO cluster sample are the Cold + Hot dark matter model 
and a low-density CDM model with a non-zero cosmological constant: all the 
other models are ruled out at a high confidence level. The pdfs of all models 
are well approximated by a lognormal distribution, even when the Gaussian 
smoothing radius is as large as 40 h~^ Mpc, consistent with the findings for 
Abell/ACO clusters. The low-order moments of all the pdfs are found to obey 
a variance-skewness relation for the form 7 Ri S'sir^, with S3 ~ 1.9, indepen- 
dent of the primordial spectrum shape and consistent with the observational 
results. Comparing the variance of the smoothed cluster density field to that of 
the dark matter distribution, we find that the linear biasing parameter for the 
simulated clusters is nearly constant over a large range of scales, but its value 
depends on the model. This suggests that it is probably a reliable procedure to 
use linear biasing to infer the dark matter power-spectrum from observational 
cluster samples. We also note that the abundances of clusters predicted using 
the Press-Schechter theory provide strong constraints on these models: only the 
mixed Cold + Hot Dark Matter model and the Iow-Hq CDM model appear to 
produce the correct number-density of clusters. Taking this constraint together 
with the cluster clustering statistics leads one to conclude that the best of our 
models is the Cold + Hot dark matter scenario. 



The study of the distribution of matter on the largest 
scales provides important constraints on models of cos- 
mic structure formation. If the gravitational instability 



picture is correct, the expected displacements of clus- 
ters of galaxies away from their primordial positions 
are much smaller than the typical separation of these 
objects. In principle, therefore, clusters of galaxies can 
yield clues about the primordial spectrum of perturba- 



2 S.Borgani, M.Plionis, P.Coles and L.Moscardini 



tions that gave rise to them. This is the reason why so 
much effort has been devoted to compiling deep clus- 
ter surveys, starting with the the pioneering work of 
Abell (1958), Zwicky et al. (1968) and Abell, Corwin 
&; Olowin (1989), and leading up to extended redshift 
surveys both in the optical (e.g. Postman, Huchra &; 
Geller 1992; Dalton et al. 1994; Collins et al. 1994, and 
references therein) and in the X-ray (e.g. Nichol, Briel 
&; Henry 1994; Romer et al. 1994) regions of the spec- 
trum. 

Accompanying the observational challenge of ac- 
quiring extended cluster redshift surveys, a great deal of 
effort has also been directed towards the provision of re- 
liable statistical characterizations of the cluster distribu- 
tion. It has been established that the cluster two-point 
correlation function is well modelled by a power-law. 



^{r) = {r/roY 



(1) 



Although the slope, 7 ~ 1.8, turns out to be quite simi- 
lar to that of galaxies, the correlation length. To, is much 
larger (cf. Bahcall &; Soneira 1983; Klypin &; Kopylov 
1983; Bahcall 1988). Different determinations, based on 
different cluster samples, indicate values in the range 
7-0 = 13-25 /i"^Mpc* (e.g., Nichol et al. 1994, and refer- 
ences therein), while more recently the reliability of the 
power-law model for ^(r) has also been questioned by 
different authors (cf. Olivier et al. 1993). 

In order to compare the observational cluster data 
sets with different cosmological models several authors 
have resorted to large N-body simulations which were 
designed to sample the length-scales relevant to the 
cluster distribution (e.g. White et al. 1987; Bahcall &; 
Cen 1992; Croft fc Efstathiou 1994). The problem with 
this kind of approach is that large N-body simulations 
are very expensive from a computational point of view. 
Therefore, one is usually forced to consider only a lim- 
ited number of models, with a small number of inde- 
pendent realizations for each model. In this respect, 
analytical approaches, based either on Eulerian linear 
theory (e.g. Bardeen et al. 1986; Coles 1989; Lumsden, 
Heavens &; Peacock 1989; Borgani 1990; Holtzman &; 
Primack 1993) or on the Zel'dovich approximation (e.g. 
Doroshkevich &; Shandarin 1978; Mann, Heavens &; Pea- 
cock 1993), are in general preferred to numerical simula- 
tions. Nevertheless, they are of limited utility since one 
is often obliged to resort to oversimplifying assumptions 
about the nature of galaxy clusters. Furthermore, statis- 
tics which go beyond the two-point correlation function 
and its Fourier transform, the power-spectrum, are hard 
to handle. Finally, it is not clear how shot-noise effects 
and/or observational biases (e.g., redshift-space distor- 
tions, selection functions, non-trivial sample geometry) 
can be realistically modelled in order to allow a consis- 
tent comparison with real data sets. 



In a previous paper (Borgani, Coles &; Moscardini 
1994, hereafter Paper I) we used the Zel'dovich approx- 
imation (ZA hereafter) to generate cluster simulations 
which were accurate when compared to N-body simu- 
lations, and, at the same time, computationally cheap 
so as to enable us to produce many realizations of sev- 
eral dark matter (DM) models (see also Blumenthal, 
Dekel &; Primack 1988). Recently, Sathyaprakash et al. 
(1994) compared the ZA and several approximations of 
non-linear gravitational clustering to direct N-body re- 
sults. They showed that, although the ZA fails to follow 
small-scale clustering in the multistream region, it is 
nevertheless able to account for non-local effects due to 
long wavelength modes in the density fluctuation spec- 
trum. For this reason, the ZA is not expected to provide 
a correct description of the internal structure and mass 
distribution of non-linear structures like galaxy clusters, 
but it is very accurate in locating them in the correct 
positions and thus reliably describes their spatial distri- 
bution. In this respect, it is not necessary to employ the 
full power and sophistication of modern N-body meth- 
ods to investigate cluster clustering on large scales. On 
the other hand, accounting for the details of non-linear 
gravitational clustering becomes crucial if one is inter- 
ested in small-scale aspects of cluster properties. 

In this paper we present cluster simulations, which 
are based on an implementation of the ZA, that has 
been substantially improved with respect to that used 
in Paper I (see also Plionis et al. 1994, hereafter Paper 
II). We will show that these simulations provide as ac- 
curate a picture of the cluster distribution as can be ob- 
tained using N-body experiments. With no significant 
computation cost we generate 50 realizations for each of 
six different initial power-spectra we will consider and 
analyze the cluster simulations using two different sta- 
tistical methods, described further below. 

The availability of such a large set of simulations 
represents an extremely powerful tool to put constraints 
on DM models through a detailed comparison of the 
statistical properties of the real and simulated cluster 
distributions. However, a reliable quantitative measure 
of cluster clustering is not easy to find. Clusters are 
rather rare objects with typical mean separation of sev- 
eral tens of Mpcs; while bright galaxies have a mean sep- 
aration comparable to their correlation length, clusters 
have a mean separation which is twice the corresponding 
To value. For this reason, shot-noise effects become im- 
portant on small scales (^ 10 h~ Mpc) while, on larger 
scales (^ 40 h~ Mpc), a low signal-to-noise ratio is ex- 
pected because the clustering is weak. Robust statistical 
estimators, which are able to provide reliable measures 
over a large range of scales, are required to describe the 
cluster distribution properly and to allow an effective 
comparison with model predictions. 



h is the Hubble constant in units of 100 km s Mpc . 



As a first statistical test, we will use in this paper 
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the quantity Js^R), which is defined through the inte- 
gral of ${r): 

MR) = i / ii^)^"<ir. (2) 

The advantage of using J3{R) over ^(r) lies on the fact 
that in a sparse distribution of objects, an integral quan- 
tity such as that defined by eq. (2), should be less sus- 
ceptible to statistical noise than a differential quantity, 
such as ${r). 

An alternative method, which is becoming increas- 
ingly popular, is the study of the probability density 
function (pdf) itself. Usually one attempts to obtain a 
continuous density field by smoothing the discrete dis- 
tribution of objects with some window function (a top- 
hat or a Gaussian one are the most commonly used). 
The smoothing procedure itself reduces significantly the 
shot-noise, which could dominate the discrete distribu- 
tion (Gaztanaga &; Yokoyama 1993). Then one can de- 
fine the pdf f{g), where g = p/{p), and derive its mo- 
ments, defined by: 



in 



-1 



S"f{S)dS, 



(3) 



where 6 = g — 1- 

The pdf and moments of different galaxy samples 
have been estimated by various authors (e.g. Saun- 
ders et al. 1991; Bouchet et al. 1993; Gaztanaga and 
Yokoyama 1993; Sheth, Mo &; Saslaw 1994; Gaztanaga 
1994). Kofman et al. (1994) have compared the pdf de- 
rived from CDM N-body simulations with that of the 
IRAS sample and the one recovered using the POTENT 
procedure with n = 1 (see also Lahav et al. 1993). Their 
main conclusion is that, if galaxies trace the mass, the 
observed pdf is consistent with Gaussian initial condi- 
tions. 

Plionis fc Valdarnini (1994, hereafter PV94) stud- 
ied the pdf (and its moments) of the 3-D Abell/ACO 
smoothed cluster distribution and compared them with 
static simulations, based on a Gaussian fluctuation spec- 
trum, which reproduced the two- and three-point clus- 
ter correlation functions as well as the observed selec- 
tion effects. They found that the real and simulated 
cluster pdf is well approximated by a lognormal dis- 
tribution. Cappi &; Maurogordato (1994) have realized 
a study of the higher-order moments for the discrete 
Abell/ACO cluster distribution, for both projected and 
redshift samples, while Kolatt, Dekel &; Primack (1994) 
estimated the pdf for real cluster samples as well as 
for cluster N-body simulations based both on Gaussian 
and non-Gaussian initial CDM fluctuations. They con- 
cluded that no evidence of non-Gaussian initial condi- 
tions are imprinted into the shape of the cluster pdf. 
In Paper II we compared the variance and the skewness 
of the smoothed Abell/ACO cluster pdf with those ob- 
tained from the ZA simulations of a list of DM models. 

In this paper we will derive the pdf statistics of our 



cluster simulations. Using the same analysis procedure 
as that used for the Abell/ACO cluster sample (PV94), 
we will be able to put stringent constraints on the mod- 
els we consider. 

The layout of this paper is as follows. In Section 2 
we describe the cluster simulations, i.e., how to optimize 
the ZA, the method of cluster identification and the con- 
sidered models for the power-spectrum. In view of the 
unsuitability of our simulation method for studying the 
small-scale structure of clusters, we also use the Press &; 
Schechter (1974; PS hereafter) method to compare the 
cluster abundances predicted by DM models with the 
available observational data. This analysis provides an 
independent constraint on the models we are consider- 
ing. In Section 3 we describe the Abell/ACO sample we 
use. In Section 4 we present the analysis of the discrete 
cluster distribution using the Js integral, while in Sec- 
tion 5 we present the pdf and moment analysis of the 
smoothed cluster distribution. In Section 6 we discuss 
our results and state our main conclusions. 



2 THE SIMULATIONS 

2.1 The Zel'dovich approach 

The Zel'dovich approximation (Zel'dovich 1970; Shan- 
darin &; Zel'dovich 1989) is based on the assumption of 
laminar flow for the motion of a self-gravitating non- 
relativistic coUisionless fluid. Let q be the initial (Eule- 
rian) position of a fluid element and r(q, t) = a(i) x(q, t) 
the final position at the time t, which is related to the co- 
moving Lagrangian coordinate x(q, t) through the cos- 
mic expansion factor a(i). The ZA amounts to assume 
the expression 



r(q,i) = a{t) [q + 6(i) Vq^(q) 



(4) 



for the Eulerian-to-Lagrangian coordinate mapping. In 
eq.(4) 6(i) is the growing mode for the evolution of lin- 
ear density perturbations and ■^(q) is the gravitational 
potential, which is related to the initial density fluctu- 
ation field, (5(q), through the Poisson equation 



(5) 



a(i) 

As a result of the factorization of the t- and q- de- 
pendence in the displacement term of eq.(4), the fluid 
particles move under the ZA along straight lines, with 
comoving peculiar velocity 



v(q,i) = i(q,i) = 6(i)Vq^(q). 



(6) 



Therefore, gravity determines the initial kick to the fluid 
particles through eqs.(5) and (6), and afterwards they 
do not feel any tidal interactions. Particles fall inside 
gravitational wells to form structures, which however 
quickly evaporate. In this sense, the ZA gives a good 
description of gravitational dynamics as far as particle 
trajectories do not intersect with each other, while its 
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validity breaks down when shell-crossing occurs, and 
local gravity dominates. 

Several prescriptions have been suggested to over- 
come the shortcomings of the ZA, such as adding a 
small viscous term to the equation of motion for the 
fluid, or by going to higher-orders in Lagrangian pertur- 
bative theory (e.g. Sahni &; Coles 1994, and references 
therein). As a further possibility. Coles, Melott &; Shan- 
darin (1993) have shown that filtering out the small- 
scale wavelength modes in the linear power-spectrum 
reduces the amount of shell-crossing, thus improving 
the performance of the ZA. Melott, Pellman &; Shan- 
darin (1993) claimed that an optimal filtering procedure 
is obtained by convolving the linear power-spectrum 
with the Gaussian filter 



WaikRf) = e 



-(kRjfl2 



(7) 



The problem then arises of suitably choosing the fil- 
tering radius Rj, in order to suppress shell-crossing as 
much as possible, without however preventing genuine 
clustering to build up. In Paper I we chose Rj so that 
the expected mass within a Gaussian window of that 
radius were of the same order (~ 10 Mq) of the mass 
for a rich galaxy cluster. The disadvantage of this ap- 
proach is that it does not rely on any objective criterion 
to optimize the ZA and treats in the same fashion dif- 
ferent fluctuation spectra, which should produce a dif- 
ferent amount of shell-crossing. The resulting filtering 
radius, Rj = 50^ h~^M.pc, is generally larger than 
the optimal ones, which we use in the present paper, 
thus causing an excessive removal of clustering. 

Kofman et al. (1994) derived the analytical expres- 
sion for the average number of streams at each Eulerian 
point, JVs, as a function of the r.m.s. fluctuation level 
of the initial Gaussian density field. In Figure 1 we plot 
JVs as a function of the filtering scale Rf for the six 
different power-spectra that we will consider (see next 
subsection), evaluated according to eq.(7) of Kofman et 
al. (1994). As a general criterion, we decided to choose 
Rf for each model so that JVs = 1.1. We found this to 
be a reasonable compromise between smaller JVs values, 
giving rapidly increasing Rf and high suppression of 
clustering, and larger JVs, at which the ZA progressively 
breaks down. The resulting r.m.s. fluctuation value cor- 
responding to JVs = 1.1 is (J = 0.88. 

By adopting this implementation of the ZA, the 
main steps of our cluster simulations are the following: 

(a) Convolve the linear power-spectrum with the Gaus- 
sian window of eq. (7), and Rf chosen as previously 
described. 

(b) Generate a random-phase realization of the den- 
sity field on 128 grid points for a cubic box of 
L = 320/i"^Mpc aside. 

(c) Move 128^ particles having initial Lagrangian posi- 
tion on the grid, according to the ZA. Each particle 
carries a mass of 4.4 x 10^^ h~^ QoMq. 



(d) Reassign the density and the velocity field on the grid 
through a TSC interpolation scheme (e.g. Hockney &; 
Eastwood 1981) for the mass and the moment carried 
by each particle. 

(e) Select clusters as local density maxima on the grid ac- 
cording to the following prescription. If dd is the aver- 
age cluster separation, then we select Nd = {L/dd) 
clusters as the Nd highest density peaks. In the fol- 
lowing, we assume dd = 40 h~ Mpc, which is appro- 
priate for the combined Abell/ACO cluster sample 
to which we will compare our simulation results (see 
Section 3). Therefore, we will analyze a distribution 
of 512 clusters in each simulation box, with periodic 
boundary conditions. 

2.2 Dark matter models 

We run simulations for six different models of the ini- 
tial fluctuation spectrum. For each model, we generate 
50 random realizations, so as to reliably estimate the 
effect of cosmic variance. All the models, except the 
open CDM one, are normalized to be consistent with 
the COBE measured quadrupole of CMB temperature 
anisotropy (Bennett et al. 1994). 

The models we have considered are the following. 

(1) The standard CDM model (SCDM), with as = 1 
for the r.m.s. fluctuation amplitude within a top-hat 
sphere of 8 h~ Mpc. 

(2) A tilted CDM model (TCDM), with n = 0.7 for 
the primordial spectral index. Tilting the primordial 
spectral shape from the scale-free one has been sug- 
gested in order to improve the CDM description of 
the large-scale structure (e.g. Cen et al. 1992; Tor- 
men et al. 1993; Liddle &; Lyth 1993; Adams et al. 
1993; Moscardini et al. 1994). 

(3) A low Hubble constant CDM model (LOWH), with 
h = 0.3. Decreasing the Hubble constant has the ef- 
fect of increasing the horizon size at the equivalence 
epoch, thus pushing to larger scales the falloff of the 
spectrum to the scale-free shape. The relevance of 
this models in alleviating several cosmological prob- 
lems has been recently emphasized by Bartlett et al. 
(1994). 

(4) A Cold + Hot DM model (CHDM), with fihot = 0.3 
for the fractional density contributed by the hot par- 
ticles. For a fixed large-scale normalization, adding 
a hot component has the effect of suppressing the 
power-spectrum amplitude at small wavelengths (see, 
e.g. Klypin et al. 1993, and references therein, for 
the relevance of CHDM). Although the small-scale 
peculiar velocities are lowered to an adequate level, 
the corresponding galaxy formation time is delayed 
so that such a model is strongly constrained by 
the detection of high-redshift objects (e.g. Klypin 
et al. 1994; Ma &; Bertschinger 1994, and references 
therein). 
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Table 1. The models. Column 2: the density parameter Hq; 
Column 3: the cosmological constant term n^l Column 4: 
the density parameter of the hot component n^ot I Column 
5: the primordial spectral index n; Column 6: the Hubble pa- 
rameter h; Column 7: the linear r.m.s. fluctuation amplitude 
at S/i-^Mpc <T8. 



Model On 



Ha 



fihot 



Rf 



SCDM 


1.0 


0.0 


0.0 


1.0 


0.5 


1.0 


4.4 


TCDM 


1.0 


0.0 


0.0 


0.7 


0.5 


2.0 


1.6 


LOWH 


1.0 


0.0 


0.0 


1.0 


0.3 


1.6 


2.4 


CHDM 


1.0 


0.0 


0.3 


1.0 


0.5 


1.5 


2.2 


OCDM 


0.2 


0.0 


0.0 


1.0 


1.0 


1.0 


4.5 


ACDM 


0.2 


0.8 


0.0 


1.0 


1.0 


1.3 


3.3 



(5) An open CDM model (OCDM), with n = 0.2. Lower- 
ing the density parameter has the effect of increasing 
the horizon size at the matter/radiation equivalence 
time, so as to shift the peak of the power-spectrum 
at larger scales. 

(6) A spatially flat, low-density CDM model (ACDM), 
with no = 0.2, Qa = 0.8 for the cosmological con- 
stant term. While neglecting the spatial curvature, as 
required by inflation, this model has been shown to 
safely describe the large-scale distribution of galaxies 
(Baugh fc Efstathiou 1993; Peacock fc Dodds 1994) 
and galaxy clusters (Bahcall &; Cen 1992; Croft &; 
Efstathiou 1994). 

The transfer functions for the above models have 
been taken from Holtzman (1989), except that of 
LOWH, which is taken from Bond fc Efstathiou (1984), 
with suitably chosen shape parameter T = Qoh = 0.3. 
All the model parameters are listed in Table 1. 

Each power-spectrum is suitably smoothed on the 
scale Rf according to the prescription described in Sec- 
tion 2.1. In Figure 1 we plot the average stream number 
per Eulerian point, JVs, as a function of the filtering ra- 
dius for all the above models. The intersection of the 
Ns = 1.1 line with each curve indicates the smooth- 
ing scale adopted for the corresponding spectrum (see 
also Table 1). Note that the larger small-scale power 
for SCDM and OCDM requires a stronger filtering to 
suppress shell-crossing. 



2.3 Cluster Abundances 

As an independent constraint on the above models, we 
have computed the expected cluster abundances, as pre- 
dicted by the standard PS formalism. If structures are 
identified through a filter W on a scale R, so as to have 
a mass M = fpR^ {p is the average matter density), 
the PS formula for the number density of objects with 
mass between M and M + dM is 



n{M)dM 



where 



1 6c r v{R) 
V2^ f Jr '^(r) 



exp 



2a^R) 






v{R) 



\R) 



27rV2(ii) 



/ 



k^ P{k) 



dW'^jkR) dk 
d{kR) kR ' 



2t^ J 



k'^ P{k)W'^{kR)dk. 



(9) 



Therefore, the total abundance of objects of mass larger 
than M is 



N{> M) 



f 

J M 



n{M')dM' . 



(10) 



In the above expressions / is a "form factor", which 
depends on the shape of the filter W: f = (27r) ' for 
a Gaussian filter, and / = 47r/3 for a top-hat filter. A 
Gaussian filter will be assumed in the following analysis. 
The parameter 6c is the critical density contrast, which 
represents the threshold value for a fluctuation to turn 
into an observable object, if evolved to the present time 
by linear theory. Arguments based on a simple spherical 
collapse suggest 6c = 1.68, but the inclusion of non- 
linear effects, as well as aspherical collapse, may lead to 
a lower value of 6c- For example, Klypin &; Rhee (1994; 
KR94 hereafter) found that the cluster mass function in 
their CHDM N-body simulations is well fit by eq.(8) by 
taking 6c = 1.5. 

White, Efstathiou fc Frenk (1993) resorted to X- 
ray data for the temperatures of the gas component 
of clusters and estimated a cluster abundance of about 
4 X 10~^{h~^Mpc)~^ for masses exceeding M = 4.2 x 
10 Mq. Using observed cluster velocity dispersions, Bi- 
viano et al. (1993) obtained an abundance of about 
6 X 10~ ( h~ Mpc)~ for clusters exceeding the above 
mass limit. 

In Figure 2 we compare model predictions at dif- 
ferent 6c values to the above observational estimates. 
Note that realistic uncertainties on cluster abundances 
are probably larger than the difference between the two 
above values. They should include variations in the av- 
erage cluster number density between different samples, 
biases toward high mass for observations of cluster ve- 
locity dispersions, uncertainties in the model used to 
relate gas temperature and cluster mass, etc. Although 
taking in mind such warnings, we note from Figure 2 
that both low-density models are ruled out, with ACDM 
producing with (5c = 1.5 more than one order of mag- 
nitude less clusters than observed. This agrees with the 
suggestion of White et al. (1993), that a higher normal- 
ization ((78 ~ 1.4) is required for these models to pro- 
duce a correct number of clusters. In the Qo = 1 mod- 
els, it turns out that the resulting abundances depend 
mostly on the as normalization value and not on the 
shape of the spectrum. This is not surprising: the Gaus- 
sian smoothing scale, R ~ 4.2 h~^M.pc, which encom- 
passes a mass of 4 x 10~ [h~ Mpc)~ , is equivalent to a 
top-hat sphere of about 7 h~ Mpc, which is rather close 
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to the normalization scale. As a result, SCDM turns out 
to produce too many clusters for any reasonable value of 
6c- On the other hand, the low normalization of TCDM 
turns into a severe underproduction of clusters, even at 
the smallest 6c values. The only two models which gen- 
erate cluster abundances in agreement with the results 
by White et al. (1993) and Biviano et al. (1993) for a 
reasonable choice of 6c are LOWH and CHDM. 

It is wise to sound a note of caution about the 
strength of the constraints emerging from the Press- 
Schechter analysis, for a number of reasons. First, there 
is some evidence of a discrepancy between the mass pro- 
files of clusters inferred from X-ray data (e.g. Edge &; 
Stewart 1991) and from gravitational lensing considera- 
tions (e.g. Kaiser et al. 1994), so it is not clear whether 
the cluster mass function inferred from the X-ray data 
is correct. Likewise, there is a possibility that the distri- 
bution of cluster peculiar velocities may be affected by 
subclustering. One should also mention that the appli- 
cability of the Press-Schechter method is itself open to 
some doubt. Although, as we mentioned above, it ap- 
pears to perform well for the CHDM model when com- 
pared with JV-body simulations, its accuracy is yet to 
be verified for the other models. In the case of open 
CDM, where a much higher fraction of the cluster mass 
is baryonic than in the other models, one might imagine 
this formalism to be particularly suspect. We therefore 
take the constraints emerging from this analysis to be 
indicative but not watertight. 



2.4 Reliability of the ZA 

Before entering into the presentation of our analysis, 
we want to stress once more the reliability of the ZA for 
simulating the large-scale distribution of galaxy clus- 
ters. In Figure 3 we plot the projected particle distribu- 
tion within a slice 10 h~ Mpc thick, superimposing the 
cluster distribution, for the SCDM and OCDM models. 
In order to better show how the identified clusters trace 
the underlying density field, we used simulations within 
a 640 h~ Mpc box^. The same initial phases have been 
used in both models, so that the structures present in 
the two slices can be directly compared. As expected, 
clusters are preferentially located at the knots corre- 
sponding to the intersections between filaments, while 
they avoid long filaments and flattened pancakes. As 
for the difference between the two models, it is appar- 
ent that the cluster distribution in the OCDM generates 
longer cluster filaments, surrounding larger underdense 
regions than in the corresponding SCDM model. Fur- 
thermore, in the OCDM simulation there are regions of 
size ~ 150 h~ Mpc which are completely devoid of clus- 



' The analysis of these simulations will be presented in forth- 
coming papers. 



ters, while in the SCDM the clusters appear to be more 
space-filling. 

In Figure 4 we compare the two-point cluster cor- 
relation function for our CHDM simulations (open dots) 
to that obtained by KR94 by evolving the same spec- 
trum with a Particle-Mesh (PM) N-body code. Our 
results refer to the average taken over 50 random real- 
izations. Error bars are estimated as the r.m.s. scatter 
over this ensemble. The KR94 results are obtained as 
an average over 2 realizations and error bars are quasi- 
Poissonian sampling uncertainties. By comparing this 
plot with Fig. 2 of Paper I, it is apparent that we have 
improved our implementation of the ZA by increasing 
the resolution and by optimizing the power-spectrum 
filtering. The agreement between ZA and N-body re- 
sults is really remarkable and extends down to quite 
small scales (~ 7 h~ Mpc), where shell-crossing should 
already play some role. This result further confirms the 
reliability of the ZA to follow correctly the mildly non- 
linear clustering regime and ensures that our simula- 
tions provide a fair representation of the cluster distri- 
bution. 



3 THE CLUSTER SAMPLE 

We use the combined Abell/ACO R > cluster sam- 
ple, as defined in Plionis &; Valdarnini (1991) [here- 
after PV91] and analysed in Plionis, Valdarnini &; Jing 
(1992) [hereafter PVJ] and Plionis fc Valdarnini (1994) 
[hereafter PV94]. The northern sample, with dec> —17° 
(Abell), is defined by those clusters that have measured 
redshift z^ 0.1, while the southern sample (ACO; Abell, 
Corwin &; Olowin 1989), with dec< —17°, is defined by 
those clusters with mio < 16.4 (note that with this def- 
inition and due to the availability of many new cluster 
redshifts only 7 ACO clusters have mio estimated red- 
shifts from the mio — z relation derived in PV91). Both 
samples are limited in Galactic latitude by |6| > 30°. 
The redshifts have been taken from a number of stud- 
ies, the references of which can be found in PV94. The 
total number of clusters in our samples is 357 and 157, 
for Abell and ACO respectively. 

To take into account the effect of Galactic absorp- 
tion, we assume the usual cosecant law: 



P(|6|) = dex [a(l -csc|6|)] 



(11) 



with a Si 0.3 for the Abell sample (Bahcall &; Soneira 
1983; Postman et al. 1989) and a w 0.2 for the ACO 
sample (Batuski et al. 1989). The cluster-redshift selec- 
tion function, P{z), is determined in the usual way (cf. 
Postman et al. 1989; PVJ; PV94), by fitting the clus- 
ter density, as a function of z (see the above reference 
for details). Cluster distances are estimated using the 
standard relation: 

^ = H./il+z) [^- + (^ - ^°)(^ - V^^^^^)] (12) 
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with Ho = 100 h km sec"^ Mpc"^ and go = flo/2. 
Strictly speaking, eq.(12) holds only for vanishing cos- 
mological constant. Therefore, for a consistent compari- 
son with the simulation models, we should use different 
R-z relations for the Abell/ACO analysis. However, we 
verified that final results are essentially independent of 
the choice of the (A, Qo) parameters used in the simu- 
lations. For this reason, in the following we will present 
results for real data only based on assuming eq.(12) with 
Qo = 0.2. 

PVJ and others have estimated the two-point cor- 
relation function for the Abell and AGO cluster samples 
and found that the slope of the two-point function has 
a value ~ 1.8 ± 0.2 for both of them. Their amplitudes, 
however, are slightly different, with Abell clusters hav- 
ing To ~ 18±4 h~ Mpc (bootstrap errors) out to $ 50 
h~^ Mpc while the AGO clusters have t-q ~ 22 ± 10 
h~ Mpc but only out to ~ 30 h~ Mpc; this fact could 
be attributed to the relatively small solid angle covered 
by the AGO and to the consequent undersampling of 
large wavelengths. Furthermore, based on the analysis 
of Jing, Plionis &; Valdarnini (1992), we believe that 
the cluster correlations are not significantly affected by 
contamination effects (cf. Sutherland 1988). 

PVJ and PV94 found that the Abell and AGO clus- 
ter number densities, out to their limit of completeness, 
are ~ 1.4x10"'^ /i^ Mpc"^ and ~ 2.1x10"'^ /i^ Mpc"^ 
corresponding to mean separations dd ~ 41 h~ Mpc 
and del ~ 36 h~ Mpc, respectively. The higher space- 
density of AGO clusters is partly due to the huge Shap- 
ley concentration (Shapley 1930), but a significant part 
is also due to systematic density differences between 
the Abell and AGO cluster samples, as a function of 
z, which has been noted in a number of studies (cf. 
PV91 and references therein) and which could be at- 
tributed to the high sensitivity of the Illa-J emulsion 
plates. In PV94, this effect was taken into account by 
normalizing the densities of the two samples using a ra- 
dial matching function, W[R), which is defined as the 
ratio between the average densities for Abell and AGO 
clusters at equal volume shells. 

In the following, we compare results based on 
the Abell/ AGO sample with those derived from our 
simulated cluster populations, selected so that dd = 
40 h~ Mpc. Variations in dd of the order of the Abell- 
AGO difference, does not significantly affect the result- 
ing statistical properties. Finally, following PV94 we re- 
strict our analysis of the real cluster sample within a 
maximum distance of Rmax = 240 h~ Mpc, in order 
to minimize the uncertainties due to the approximate 
character of the redshift selection function, P{z), and of 
the radial maic/img function, W[R), especially at large 
distances. 



4 STATISTICS OF THE DISCRETE 
CLUSTER DISTRIBUTION 

Our first statistical test for comparing the real data and 
the simulations, involves the evaluation of the quantity 
J3{R), defined by eq. (2). It is straightforward from the 
definition of this quantity to construct the estimator 



MR) 



R^ /N, 
3 



(f^-). 



(13) 



where JV„{, is the average number of cluster neighbours 
within a distance R from a cluster, while N is the ex- 
pected number of neighbours for a random cluster distri- 
bution (estimated at the positions of the real clusters). 
Therefore, J3{R) oa R 'as long as ${r) oc r '. 

It has been argued (cf. KR94) that the scale at 
which the power-law shape of ^(r) breaks and firstly 
crosses zero, is a potentially powerful test for cosmo- 
logical models. However, since such a scale corresponds 
by definition to the weak clustering regime, its detec- 
tion can be heavily affected by statistical noise; for an 
explicit demonstration of this, see Paper I. In this re- 
spect, the analysis of Js should have the advantage of 
being more stable and suffering less from observational 
biases. 

In Figure 5 we plot J3{R) for the real data (filled 
circles) and simulations (open circles). We estimated N 
for the real data by averaging over 100 random samples, 
having the same selection criteria (boundaries, galac- 
tic extinction function, redshift selection and system- 
atic Abell/ AGO differences) as the real one. Error bars 
for the simulated samples are la scatter over the en- 
semble of 50 realizations. In Table 2 we report values 
of Ja for data and simulations at three different scales. 
The quoted uncertainties for real cluster analysis are la 
scatter estimated over an ensemble of 100 bootstrap re- 
samplings. Such errors are not plotted in Figure 5. In 
fact, since we are asking which is the probability that 
a given model generates a result like that of the ob- 
served cluster distribution, its 'success' is just measured 
by the distance of the real data point from the cos- 
mic r.m.s. error bars. This should be taken into account 
when judging to which confidence level a model has to 
be accepted or rejected. In Paper II we verified that in- 
termediate scales of few tens of Mpcs are best suited to 
constrain DM models, when using the cluster distribu- 
tion, smaller and larger scales being affected by shot- 
noise and low signal-to-noise ratio, respectively. In the 
present analysis we do not consider scales much smaller 
than 20 h~ Mpc as well as larger than 60 h~ Mpc. 

For the Abell/ AGO sample, J3{R) increases up to 
ii ~ 35 h~ Mpc, flattens at a scale corresponding to 
the break of the power-law shape of ${r), and even- 
tually declines at R^ 50 h~ Mpc, after which ${r) be- 
comes negative. By comparing this result with those of 
the simulations, it turns out that the only two mod- 
els which overcome this test are GHDM and AGDM, 
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Table 2. Js[R) values at different scales (in h Mpc units) 
for both simulated and real cluster distributions. 



Model 



28.9 



i « 43.2 



64.3 



SCDM 5.00 ±0.64 5.55 ± 1.34 3.86 ± 2.12 

TCDM 4.93 ±0.74 6.05 ± 1.60 5.71 ± 2.74 

LOWH 5.12 ±0.79 6.10 ± 1.60 5.16 ± 2.83 

CHDM 8.39 ±0.94 10.19 ± 1.98 8.77 ± 3.36 

OCDM 11.89 ±1.73 16.52 ± 2.80 18.66 ± 6.05 

ACDM 8.43 ±1.03 10.99 ± 2.19 11.61 ±4.27 

Abell/ACO 8.10 ±1.69 10.28 ± 2.53 6.96 ± 5.17 

although both of them seem to produce too strong clus- 
tering at the smallest scale considered and the second 
one is also marginally overclustered at the largest scale. 
All the other models are ruled out at least at a 3a level. 
The cluster distributions for the SCDM, LOWH and 
TCDM models are too weakly clustered over the whole 
scale range. Note that the SCDM has a rather flat J3{R) 
profile, according to the expectation that this model has 
a cluster two-point correlation function which declines 
rapidly beyond ~ 20 h~^M.pc. Conversely, OCDM gen- 
erates too much clustering, with J3{R) increasing up 
to R^ 60 h~ Mpc. These results confirm what we al- 
ready found in Paper II from the analysis of the vari- 
ance/skewness relation of the smoothed cluster density 
field. Note that the results for TCDM and LOWH are 
remarkably similar. This agrees with the expectation 
that, as far as the shape of the power-spectrum is con- 
cerned, a change in the Hubble parameter h is roughly 
equivalent to a change in the spectral index n according 
to the relation Ah = -An (cf. Lyth fc Liddle 1994). 



5 STATISTICS OF THE SMOOTHED 
CLUSTER DISTRIBUTION 

In order to facilitate the comparison of our results with 
those obtained by PV94 from the combined Abell/ACO 
cluster sample, we followed basically the same procedure 
as they did, and which we briefly describe below. We 
obtain a continuous cluster density field by smoothing 
the cluster distribution on a grid, with grid-cell width of 
20 h~ Mpc (16 grid-points), using a Gaussian kernel: 

W(|xi-xg|)= (27riiL)"''' exp 



^gl 



2R-i 



.(14) 



The smoothed cluster density, at the grid-cell positions 
Xg, is then: 



P(xg) 



E,p(x.)yv(|x.-xg|) 



(15) 



/W(|x-xg|)d^^ ' 

where the sum is over the distribution of clusters with 
positions x;. In order to study the cluster density field 



at different smoothing scales, we use three radii for the 
Gaussian kernel: Rsm = 20, 30 and 40 h~ Mpc with 
|xj- — Xgl < SRsm- Therefore the integral in the denomi- 
nator of eq. (15) has a value smaller than unity (~ 0.97). 

5.1 The probability density function 

As a first test for the smoothed cluster density field, we 
work out the probability density function, f{g), which 
represents a low-order (one-point) statistics. We then 
compare the pdf of each set of cluster simulations with 
the observed Abell/ACO pdf, derived by PV94, as well 
as with the following theoretical models. 



(a) The Gaussian distribution given by 



/(e) 



V^ 



: exp 



(g-lf 



2cr2 



(16) 



where a is the standard deviation of g (= p/{p)). If 
f{g) is a Gaussian then it should be defined in an 
infinite interval, which implies that f[g < 0) /: 0. 
Since, however g > hy definition, f{g) is expected 
to be well approximated by a Gaussian only in the 
limit (J ^ 0. In this case the skewness, 7 (= {(5^)), 
vanishes. Even in the case of an initial Gaussian den- 
sity field, the gravitational evolution acts in such a 
way as to increase the variance a , and thus, due to 
the constraint g > 0, f{g) has to become positively 
skewed. For as long as the variance a = {6 ) is small, 
the deviation of the pdf shape from a Gaussian is well 
approximated by the Edgeworth expansion (Colombi 
1994). 
(b) The lognormal distribution given by 

1 r (Ine-MiT' 



m 



V^ 



: exp 



2al 



(17) 



where g is obtained through an exponential trans- 
formation of a Gaussian random variable % as g = 
exp(x). In eq. (17), /il and aL are the mean and 
standard deviation of In g respectively. It has been ar- 
gued that this distribution describes the distribution 
of density perturbations resulting from Gaussian ini- 
tial conditions in the weakly non-linear regime Coles 
&; Jones 1991). Bernardeau &; Kofman (1994) have 
shown that the lognormal distribution is not really 
a natural consequence of mildly non-linear gravita- 
tional evolution, but a very convenient fit only in some 
portion of the [a, 7i)-plane (i.e. cr <C 1 and spectral in- 
dex n Si —1). It has nevertheless been found to give an 
extremely good fit to the CDM density and the IRAS 
galaxy pdf in the weakly-linear regime (Kofman et al. 
1994), as well as to the observed Abell/ACO cluster 
distribution (PV94). 
(c) The pdf following from the Zel'dovich approximation 
(Kofman et al. 1994): 

fU) = JJil 

■'^^' 47rJV3eV 
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poo 



dse-^'-''>"'"'Ul+e- 



2 /i^S 



X lei' + e 

/3„(s) = -y5s{l/2 +cos[2/3 (ti- l)7r 
+ 1/3 arccos (54/es^ - l)]}, 



3/20- 






(18) 



where a is the rms amplitude of density fluctuations 
and Ns is the average stream number per Eulerian 
point. 

Note that discreteness effects could be important 
since they affect the shape of the pdf and the estimation 
of its moments; especially at small Rsm, when the num- 
ber of clusters in the Gaussian sphere is small, and/or 
when the smoothing fails to create a continuous density 
field due to discreteness (in our case this is apparent in 
the Rsm = 20 h~^ Mpc case for g < 0.8). In the case of 
a Poisson sampling of an underlying continuous density 
field, the shot-noise contributions to the moments can 
be easily estimated and corrected for (cf. Peebles 1980). 
However, the cluster distribution can hardly be consid- 
ered a Poissonian sampling of the underlying (galaxy) 
distribution, since clusters are expected to form at high 
density peaks. Therefore, the Poissonian shot-noise cor- 
rection could not give a reasonable description of dis- 
creteness effects (Coles &; Frenk 1991; Borgani et al. 
1994). Moreover, Gaztanaga &; Yokoyama (1993) have 
shown that the smoothing process itself considerably 
suppresses these shot-noise effects. For these reasons 
PV94 did not use any shot-noise corrections. To make 
a consistent comparison of our models with the data, 
we also did not include such corrections in our analysis. 
Since all the model cluster distributions have the same 
mean number density and we treat them similarly, the 
possible effects of shot-noise are accounted for in the 
same way in both the data and the simulations: we are 
therefore comparing like with like. 

In Table 3 we present the results of the compari- 
son between the simulation pdf, the PV94 Abell/ACO 
cluster pdf, the lognormal as well as the Gaussian dis- 
tributions, using a % -test defined as: 



X 



E 



fris) - /.■''°°^(g) 



(19) 



where the weights e? correspond to cosmic variance. 
Note that we derived the simulation pdf in redshift 
space so that a consistent comparison with the PV94 
results can be made, while the comparison with theo- 
retical models is done in real space. 

In Figure 6 we present the simulation cluster pdfs 
for the most successful and least successful models 
(i.e. CHDM and SCDM), together with the PV94 
Abell/ACO pdf at Rsm = 20 and 40 /i"^ Mpc, and in 
Figure 7 we make the comparison between these two 
cosmological models and the theoretical distributions. 
The error bars represent the scatter around the ensem- 



Table 3. x probabilities that the indicated simulation 
model pdf could have been drawn from a parent distribution 
given by the lognormal, Gaussian or the real cluster (PV94) 
pdfs. For the Rsm = 20 h~ Mpc case the comparison is 
performed for g > 0.8 (see text). Also reported in Column 6 
are the values of the reduced skewness ^'3 for real simulations 
and for real data. 



Model Rs 



SCDM 



TCDM 



LOWH 



CHDM 



OCDM 



ACDM 



30 
40 

20 
30 
40 



30 
40 

20 
30 
40 



0.83 
0.99 

0.15 
0.99 
0.99 



0.89 
0.99 

0.10 
0.99 
0.99 



0.00 
0.92 



0.75 
0.97 



0.00 4 X 10" 
0.00 0.99 
0.06 0.99 



0.00 
0.49 

0.00 
0.00 
0.00 



20 3.8 X 10"^ 0.00 



0.97 
0.97 

0.95 
0.98 
0.63 

0.68 



1.83 ±0.70 
1.60 ±1.25 

1.96 ±0.32 
2.02 ±0.80 
1.96 ±1.44 



1.87 ±0.69 
1.75 ±1.18 

1.93 ±0.30 
1.96 ±0.65 
1.81 ±1.08 

1.87 ±0.21 



30 
40 

20 
30 
40 



0.94 
0.99 

0.39 
0.33 
0.02 



0.00 7.5 X lO"-" 1.87 ±0.41 
0.00 1.1 X 10"^ 1.77 ±0.72 



0.00 0.98 
0.00 0.72 
0.70 0.22 



20 2 X lO"-* 0.00 

Abell/ACO 30 6.3 x 10"^ 0.00 

40 0.99 0.27 



1.90 ±0.26 
1.96 ±0.55 
1.95 ±0.91 

1.81 ±0.23 
1.78 ±1.30 
1.76 ±1.85 



ble mean values (cosmic variance). As in Table 3, com- 
parisons with theoretical models and real data are made 
in real space and in redshift space, respectively. There 
is an excellent agreement between the CHDM and the 
Abell/ACO cluster pdfs while in the SCDM case there 
is a clear discrepancy at small (< 0.5) and large g's. 
It is apparent that: 

(i) The Gaussian distribution does not provide a good 
fit at any Rsm < 30 h~ Mpc and for any model. 
For Rsm = 40 h~ Mpc the Gaussian fit is accept- 
able only for the SCDM, LOWH and ACDM mod- 
els, 
(ii) The Zel'dovich pdf model is inconsistent with the 
simulation results, even though the underlying dy- 
namics governing the cluster distribution are de- 
scribed by the ZA. We find that this distribution is 
ruled out at a confidence level larger than 99.99%, 
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for Rsm < 30 /i Mpc for all the simulation mod- 
els. This is the reason why we did not show re- 
sults for this model in Table 3. The SCDM, TCDM, 
CHDM and OCDM models are only consistent at 
a ~ 20% - 25% level for Rsm = 40 h'^ Mpc. One 
may argue that, since the ZA pdf is designed to 
describe the DM clustering, its failure for the clus- 
ter distribution is nothing but the consequence of 
not accounting for the mass within clusters in the 
analysis of their distribution. To check this, we re- 
peated the analysis by weighting simulated clusters 
according to their mass and found no appreciable 
differences in the pdf shapes. 

(iii) All simulation pdf 's are well approximated by a log- 
normal distribution, irrespective of their different 
power-spectra. Note that at the Rsm = 20 h~ Mpc 
case the comparison is done for g > 0.8 because, 
at lower values of g, discreteness effects introduce 
significant noise. Therefore, in contrast to the case 
of the matter distribution, the lognormal fit to the 
cluster pdf is more likely to be connected with the 
high-peak biasing description of cluster formation 
than being due to non-linear gravitational effects, 
which dominate much smaller scales (Bernardeau 
fc Kofman 1994). 

(iv) The scale which best discriminates between dif- 
ferent models and the Abell/ACO data is clearly 
Rsm = 20 h~ Mpc. At larger scales, all the models, 
except OCDM, produce acceptable fits. The only 
models that produce a pdf consistent with the P V94 
results, at all 3 smoothing radii, are the CHDM and 
the ACDM, although CHDM fares better on larger 
scales, in agreement with the J3{R) results and with 
the results of Paper II. 



5.2 Moments of the pdf 

According to eq. (3), the moments of the pdf give a 
large weight to the high density tail ((5 > 1) of the pdf. 
They are therefore expected to suffer less from shot- 
noise effects, which are smaller in the overdense parts of 
the distribution. In Paper II we presented results about 
the variance, a"^ = (<5^), and the skewness, 7 = (<5^). 

It has been argued on several grounds (e.g. Coles 
&; Frenk 1991) that the relation 

J^Ssia'Y, (20) 

with 53 nearly independent of scale, should describe 
the clustering of cosmic structures. Although at small 
scales, below a few Mpc, eq. (20) is predicted by mod- 
els of non-linear gravitational clustering (e.g. Borgani 
1994 and references therein), at the larger scales, sam- 
pled by galaxy clusters, it is expected to hold due to 
mildly non-linear evolution as well as by the bias relat- 
ing the cluster and DM distributions. The resulting S3 
values at different Rsm are reported in Table 3 for both 



the simulations and the Abell/ACO sample. We note 
that (a) the reduced skewness S3 is always independent 
of the scale with a good accuracy, and (b) it takes the 
same value S3 ~ 1.9 for all the models, within statistical 
fluctuations, and consistent with the observational data 
(PV94). 

Accordingly, we conclude that, for the cluster distri- 
bution, only the amplitude of clustering, and not its na- 
ture, depends on the initial power-spectrum. This sug- 
gests that both the lognormal pdf shape and the S3 
value observed for real data are natural consequences 
of high-peak biasing and possibly of the random-phase 
assumption of the primordial density field. 



5.3 The bias parameter of clusters 

It is interesting to attempt to derive from our simula- 
tions the cluster biasing parameter, which is defined as 
the ratio between the r.m.s. fluctuations of the cluster 
and DM distributions: 



(TDM 



(21) 



The relevance of the biasing parameter lies in the fact 
that it is a necessary ingredient if one is seeking to de- 
rive the shape of the spectrum of the initial density 
fluctuations, starting from the distribution of clusters 
of galaxies. The usual assumption is that biasing is lin- 
ear, that is bci does not depend on the scale. From their 
power-spectrum analysis, Jing &; Valdarnini (1993) and 
Peacock &; Dodds (1994) verified that at least the rel- 
ative biasing between clusters and optically as well as 
infrared selected galaxies is independent of the scale to 
a quite good accuracy. In particular. Peacock &; Dodds 
(1994) found that bci/bjRAS — 4.5, for the relative bias- 
ing between clusters and IRAS galaxies. In our simula- 
tions, we can directly verify whether the linear biasing 
paradigm for clusters is reasonable, since we also know 
the variance of the linear DM fluctuations: 



0'{Rsm)Dl\ 



IT^ J 



dkk"^ P{k)WR,^{k) 



1/2 



,(22) 



where VVii,„(A;) = exp( — A; Rsm/^) is the Fourier trans- 
form of the window function of eq. (14). 

In Figure 8 we plot the biasing parameters for the 
cluster simulations at different Rsm- The heavy-dashed 
horizontal lines delineates the la band for the observa- 
tional result, once we assume that IRAS galaxies fairly 
trace the DM distribution, bjRAS = 1.0 ± 0.2. Indepen- 
dent of the model considered, bd is fairly constant over 
the whole scale range. Figure 8 is quite different from 
the analogous plot presented in Paper I, which showed 
a decreasing trend for bd at small scales. This confirms 
how important is the increased resolution and the opti- 
mization of the ZA for the reliability of our simulations. 
From one hand, such a linearity of the biasing is a rather 
remarkable result, since both the evolution of the den- 
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sity field and the selection of clusters as high-density 
peaks represent definitely non-linear transformations of 
the initial fluctuations. From the other hand, this result 
supports the usual assumption of bd = const, used to 
infer the shape of the primordial power-spectrum from 
that of clusters. Note, however, that different models 
have rather different values of bd . The TCDM clusters 
are by far the most biased. The resulting bd — 7 could 
be made consistent with observations only by allowing 
for a IRAS biasing as low as bjRAS — 0.56. 



6 CONCLUSIONS 

In this paper we have compared results of a statistical 
analyses of our simulated cluster distributions and of 
a combined sample of Abell/ACO clusters. Our cluster 
simulations, which are based on the Zel'dovich approx- 
imation, are extremely cheap computationally; each re- 
alization takes about 6 minutes of CPU on a IIP755/125 
workstation. This has allowed us to run a large number 
of realizations (50) for each model, so as to properly es- 
timate the cosmic variance. We have also carefully ver- 
ified that such a low computational cost is not at the 
expense of reliability. Our cluster selection procedure 
identifies clusters in correspondence of knots at the in- 
tersection of filaments (see Figure 3). Comparing the 
resulting two-point correlation function for our CHDM 
clusters with that obtained by Klypin &; Rhee (1994) 
from their PM simulations of the same model, we find 
a extremely good agreement, even down to quite small 
scales (~ 7 h~ Mpc), at which one can doubt the valid- 
ity of the ZA. The reason for this remarkable success of 
the ZA in reproducing N-body results lies in its abil- 
ity to account for non-local effects when moving parti- 
cles from their initial (linear) positions to their correct 
evolved ones (cf. Pauls &; Melott 1994; Sathyaprakash 
et al. 1994). 

As a preliminary indication of the success of the 
DM models, we compared the predicted cluster abun- 
dances with available observational results using the 
Press &; Schechter (1974) approach. The only models 
that produce adequate abundances are CHDM and low- 
Ho CDM but, given the uncertainties in both theory and 
observation, we take these results to be indicative rather 
than definitive. 

We have also analyzed both the point-like cluster 
distribution using the Js statistic, and the Gaussianly- 
smoothed cluster distribution. Both analyses provide 
consistent stringent constraints on the simulated DM 
models, which confirm the results presented in Paper II. 
The only models which pass all the tests are the Cold + 
Hot DM model and the low-density CDM model with 
non-vanishing cosmological constant (ACDM). Stan- 
dard, Tilted and Low-Ho CDM versions do not account 
for the large-scale clustering of the real cluster sample. 



The Open CDM model produces a much stronger clus- 
tering than real clusters on all scales. 

We find that the shape of the probability density 
function (pdf) is not well approximated by a Gaussian 
distribution, even at the largest smoothing scale con- 
sidered. The pdf for both simulations and Abell/ACO 
clusters is always much better reproduced by the lognor- 
mal model than by the Zel'dovich prediction, despite the 
fact that the ZA governs the underlying dynamics. This 
shows that, at least for clusters, the lognormal shape of 
the pdf does not occur by chance, for a limited set of 
initial conditions, as argued to happen for the galaxy 
distribution (Bernardeau &; Kofman 1994). Instead it is 
much more likely to be related to the fact that clusters 
trace the high density peaks of the underlying matter 
field and, perhaps, to the initial random-phase assump- 
tion. Furthermore, using a % -test to compare the shape 
of the data and simulations pdf, we find that the best 
models are the CHDM and ACDM ones, in agreement 
with the Ja and the moment analysis (Paper II). 

Although the variance and skewness of the 
smoothed cluster pdf are powerful discriminators of dif- 
ferent models (see Paper II), the reduced skewness, 
S3 = T/cr , turns out to be independent of the initial 
spectrum. We always find S3 ~ 1.8-1.9, almost inde- 
pendent of the scale, and consistent with the observa- 
tional results (Plionis &; Valdarnini 1994). One could be 
tempted to conclude that such a value of S3 is naturally 
produced by the high-peak selection of clusters, prob- 
ably combined with the Gaussian nature of the initial 
fluctuations. Whether the analysis of the reduced skew- 
ness represents a test of Gaussian vs. non-Gaussian ini- 
tial conditions remains to be seen (see also Coles et al. 
1993). 

We have verified that the linear biasing prescrip- 
tion used to relate cluster and DM distributions is al- 
ways satisfied to a good precision. The resulting value 
of the biasing parameter, bd, depends on the details of 
the model, so that it could represent a further discrim- 
inator between different initial power-spectra. Taking 
bd/biRAS — 4-5 for the ratio between the biasing pa- 
rameters for clusters and IRAS galaxies (cf. Peacock &; 
Dodds 1994; Plionis 1994), and assuming that clusters 
trace the DM distribution fairly, we find that OCDM 
and LOWH models have a marginally too low and too 
high bd , respectively, while TCDM clusters are far too 
strongly biased tracers of the density field, compared to 
real clusters. 

The overall picture emerging from these studies is 
that the large-scale cluster distribution places stringent 
constraints on models of structure formation. Observed 
cluster clustering is rather well reproduced either by the 
CHDM model with fihot = 0.3 or by a low density CDM 
model no = 0.2, with non-vanishing cosmological con- 
stant term, Qa = 0.8. If combined with results on the 
cluster abundances, this suggests that the only surviving 
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model is CHDM, with ACDM producing more than one 
order of magnitude fewer clusters than observed. Fur- 
thermore, since ACDM behaves much like OCDM from 
the point of view of the large-scale velocity fields, our 
hope is that we should be able to further discriminate 
between CHDM and ACDM by analyzing the dipole 
structure of the cluster distribution (Tini Brunozzi et 
al. 1994) and cluster peculiar velocities. 

One can also ask whether reasonable modifications 
of the parameters in the DM models considered above 
(i.e. the Cold + Hot DM mixture, the no and A values, 
the primordial spectral index n, the Hubble constant, 
etc.) could lead to significant changes in the resulting 
cluster distribution. Verifying this will not be a difficult 
task, thanks to the low computational cost of our sim- 
ulations, which makes them a fiexible instrument to ex- 
plore the entire parameter space of DM models. We be- 
lieve that, in future investigations of specific DM mod- 
els, the optimal strategy would be, first of all, to run 
optimal ZA simulations in order to assess the model on 
large scales. Only after that one should decide whether 
a model is worth exploring at smaller scales by means 
of high-resolution, computationally expensive N-body 
simulations. 
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Figure Captions 

Figure 1. The average number of streams per Eule- 
rian point, JVs, as a function of the Gaussian filtering 
scale Rf for the six different models. The dotted hori- 
zontal line, JVs = 1.1, delineates the quasi single-stream 
regime. 

Figure 2. The abundances of clusters with mass M > 
4.2 X lO'^^M© for the different DM models, as predicted 
by the Press &; Schechter (1974) formalism, as a func- 
tion of the critical density contrast 6c- The horizontal 
lines are the observational results by White et al. (1993; 
dotted line) and by Biviano et al. (1993; dashed line) 
for clusters with mass larger than the above value. 
Figure 3. The cluster distribution (heavy dots) su- 
perimposed on the DM particle distribution in a slice 
10 h~^Mpc thick for a box of side 640 h~^Mpc. The up- 



per panel is for the SCDM model and the lower panel for 
the OCDM model. Initial random phases are the same 
for the two models, so that the structures in the plots 
can be compared directly. It is interesting to note that 
clusters are strongly correlated with the intersection of 
filaments of the DM distribution. 

Figure 4. Comparison between the two-point correla- 
tion function for our CHDM cluster simulations (open 
dots) and for the PM N-body results from Klypin &; 
Rhee (1994) based on the same initial spectrum (filled 
dots). Their results are obtained with PM simulations 
having 256 grid points on a box of side 200 h~ Mpc. 
Our results are the average over 50 realizations, and er- 
rors are la scatter over this ensemble. The KR94 results 
are an average over 2 realizations and the error bars are 
quasi-Poissonian estimates. 

Figure 5. The J3{R) integral as a function of the scale 
R for the simulated (open circles) and real Abell/ACO 
(filled circles) cluster distributions. Error bars are plot- 
ted only for the simulations and correspond to la scatter 
over the ensemble of 50 realizations. 
Figure 6. Comparison between the pdf 's for real (filled 
circles) and simulated (open circles) cluster distribu- 
tions. Results only for the SCDM and CHDM are shown 
at Rsm = 20 and 40 h~ Mpc. For the simulations, the 
analysis is realized in redshift space and the error bars 
correspond to cosmic r.m.s. scatter. 

Figure 7. Comparison between the pdf 's of simulated 
cluster distributions and the theoretical models. We plot 
only results for SCDM and CHDM models at Rsm = 20 
and 40 h~^M.pc. Solid, long-dashed and short-dashed 
curves correspond to the lognormal, Zel'dovich and 
Gaussian model, respectively. Error bars are cosmic 
r.m.s. scatter. 

Figure 8. Scale dependence of the cluster biasing pa- 
rameter for the six different models. For reasons of clar- 
ity, we plot the corresponding cosmic r.m.s. scatter only 
for the LOWH model. Similar uncertainties hold also 
for the other models. The heavy dashed horizontal lines 
show the range indicated by observational results (see 
text). 
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